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Abstract 

We construct small-world spring networks based on a one dimensional chain and study its static 
and quasistatic behavior with respect to external forces. Regular bonds and shortcuts are assigned 
linear springs of constant k and k', respectively. In our models, shortcuts can only stand extensions 
less than 6 C beyond which they are removed from the network. First we consider the simple cases of 
a hierarchical small- world network and a complete network. In the main part of this paper we study 
random small- world networks (RSWN) in which each pair of nodes is connected by a shortcut with 
probability p. We obtain a scaling relation for the effective stiffness of RSWN when k = k' . In this 
case the extension distribution of shortcuts is scale free with the exponent —2. There is a strong 
positive correlation between the extension of shortcuts and their betweenness. We find that the 
chemical end-to-end distance (CEED) could change either abruptly or continuously with respect 
to the external force. In the former case, the critical force is determined by the average number of 
shortcuts emanating from a node. In the latter case, the distribution of changes in CEED obeys 
power laws of the exponent —a with a < 3/2. 
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I. INTRODUCTION 



Recent studies indicate that real networks have a complex structure and function 
[H 3, Q J, si- The small- world network (SWN) introduced by Watts and Strogatz [3] 



captures some basic ingredients of real networks. There are lots of studies dealing with 
static properties of SWNs, see [3| and references therein. However, we still only know a 
little bit about dynamical features of small- world networks jf], ^SSC U|- To best of our 



knowledge there is no study on the elastic properties of SWNs. Certainly the most physical 
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ication of small-world networks is in the context of macromolecules and polymers 
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id ]. The interaction pattern of a polymer can be represented by a SWN. In 



this way, the monomers are mapped into nodes of a network and the interaction between 
two spatially close monomers is shown by a bond between the corresponding nodes. Note 
that here there are two kinds of bonds: regular and long range. The long range bonds that 
are usually weaker than the regular ones, could stand much smaller stresses and would 
be broken more easily. If we model the interactions with linear springs, we will obtain 
a small-world spring network. Of course, we should take more care on modeling linear 
polymers with SWNs Nevertheless, we expect that the study of small- world networks 
(as a toy model of polymers when thermal activities are absent) provides insights about 
more complicated behaviors of real polymers. 

Suppose that we have a chain of elastic fibers and add some other fibers to randomly 
selected pairs of nodes. Then it is interesting to know , for example, the effective stiffness of 
this object and some other quantities which are of interest in fiber-bundle models }l8. 19|. 
In this paper we are going to study the behavior of small-world spring networks when 
external force F is exerted on the end nodes. The network response contains some infor- 
mation about its internal structure. For a given F, we obtain the network stiffness and the 
extension distribution of long range bonds (shortcuts). Then we increase F quasistatically 
and define a cutoff length for shortcuts beyond which we remove them from the network. 
During this quasistatic process the chemical end-to-end distance (CEED), defined as the 
number of bonds in the shortest path connecting node 1 to N, could have nontrivial 
behavior with F. Moreover by increasing F one encounters a number of avalanches in 
which one has a change in the number of shortcuts and also in CEED. Certainly, distinct 
structures could lead to different behaviors, for instance in the distribution of changes in 



2 



CEED. These differences provide us a useful measure to classify various networks. 
Here we will show that the extension distribution is scale free in the small-world regime in 
contrast to the three-modal distribution of complete-network regime. We obtain a positive 
correlation between the betweenness Q] of a shortcut and its extension. There is also a 
threshold value F c which leads to an abrupt change in CEED. The scaling of F c with the 
size of the network depends strongly on the network structure. Depending on the elastic 
properties of springs we could also have a continuous change in CEED. 

The structure of this paper is as follows. After giving some general definitions, we in- 
troduce and study a hierarchical small- world network in section Hm In section HV1 we study 
complete networks as another simple case. In section |3 we present the results of our numer- 
ical simulations and scaling arguments for random small-world networks . Finally we give 
the conclusion remarks. 

II. GENERAL DEFINITIONS 

We take a one dimensional chain with N nodes, numbered 1 to N. This chain has N — 1 
regular bonds, with spring constant k. The network structure is completed by adding some 
shortcuts of spring constant k'. In the following we will consider the case < k' < k 
which is more reasonable in physical models. The effective elastic constant of the network 
is denoted by K. We exert force F on the end nodes of a chain and obtain the extension 
distribution of shortcuts, P(A X ), where k'A x is the force acting on a shortcut. On increasing 
F quasistatically, we define a cutoff length 6 C for shortcuts. It means that a shortcut is teared 
if its extension exceeds 5 C . Here we will take 5 C — 1 (the specific value of 5 C does not affect 
the qualitative behavior of our results). For a given network, we may have a number of 
avalanches in some particular forces. An avalanche starts by the tearing of a shortcut and 
ends when all the shortcuts have an extension smaller than 5 C . By increasing F, a number 
of avalanches may occur which we label by index a. In each avalanche, one can measure, for 
example, the force value F a , the chemical end-to-end distance R a and the effective spring 
constant of the network K a . The CEED, R, is defined as the number of bonds in the 
shortest path connecting node 1 to N. Just after the a— th avalanche, the physical end-to- 
end distance, X a , is given by F a /K a . We also define the change in a quantity such as R by 
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FIG. 1: Constructing a hierarchical small-world network. 



AR a :— R a+ i — -Ra- 
III. A HIERARCHICAL SMALL- WORLD NETWORK 

Let us first study a simple hierarchical small-world network where its behavior can be 
treated exactly. The construction of the network is depicted in Fig. [TJ In the first step, 
t — 1, we consider two nodes connected with a regular bond. In the next step, t = 2, we 
make a copy of the previous step and merge them. We also add an additional shortcut 
between the first and last nodes of the new network. If we repeat this procedure for t steps 
we obtain a hierarchical SWN of size N(t) = 2 1 - 1 + 1 with M(t) = 2 1 - 1 - 1 shortcuts. One 
can obtain the following relation for the effective stiffness of the network 

K (f)=fc . + ^ = _!_ +2 ,< (1 __i_ ) . (1) 

If we pull the end nodes of the network with force F we find that there are 2 l shortcuts of 
the extension A x = F/(K(t)2 l ). Here I e [0,t - 2] is an integer that labels the shortcuts 
according to the step they have been added to the network. This in turn results in the 
following extension distribution 

P(A X ) ~ (A x )~ 2 . (2) 

Now we start to increase F from zero. In very small F, Rq = 1. When the force reaches 
to Fx = K(t)5 c , the spring connecting node 1 to N(t) extends by S c and tears. After this 
event, R± = 2 and the number of shortcuts decreases to Mi = M(t) — 1. Moreover the 
effective network stiffness is given by K\ = Kit — l)/2. It is easy to see that a shortcut 
connecting two far nodes has a larger extension than one connecting two near nodes. Thus 
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we can summarize the behavior of interesting quantities versus a 



F a = K{t-a + l), (3) 

Ra = 2 a , 

M a = M(t) + 1 - 2 a , 

The physical end-to-end distance then reads 

X. * = B.^t±3., (4) 
K a K[t - a) 

which for t — a ^> 1 is approximately equal to R a . Thus in this limit the physical end-to- 
end distance is well given by CEED. From the above relations we obtain F a = K a -.\R a -\. 
Consider the bonds of the shortest path connecting node 1 to N. F a is the force at which 
the extensions of these bonds become S c . Note that, however, the relation between force and 
CEED is linear 



^ + 2k'(l - -' 
2* v 2 

From Eqs. © one easily finds that 



F a = ^ + 2k\l-^). (5) 



k IV 

AR a = —AM a = 2 a , AF a = — (1 - — ). (6) 

We see that for k' < k/2, F a increases versus a. Thus in this situation we have a continues 
change in R. In the continuum approximation we obtain 

P(AR) ~ (AR)~\ (7) 

On the other hand, when k' > k/2, from Eqs. (JHJ) and (jOJ we have F\ = K(t) and F\ > F2 > 
F3 . . .. It means that after tearing the first shortcut, we reach a new configuration in which 
the extension of most extended shortcut (s) is(are) greater than S c . This process continues 
until all shortcuts tear. Note that during this event which occurs at F — F\, CEED changes 
from 1 to AT — 1. Thus we have an abrupt change of CEED at F c = K(t). 



IV. COMPLETE NETWORKS 



In a complete network, each node is connected to all the other nodes. Thus the number 
of shortcuts is M = (N — 1)(N — 2)/2. We start with the simple case of k! — k with the 
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FIG. 2: Chemical end-to-end distance vs a for a complete network of size TV = 1000 with k = 1 
and k' = 0. 

effective network stiffness of K = kN/2. Moreover, P(A X ) consists of three delta peaks 
corresponding to three different kinds of shortcuts in the network: The first kind is a single 
shortcut between the end nodes and has the largest extension. The second kind of shortcuts 
which have less extensions, are those that connect the end nodes to the inner nodes. And 
the remaining shortcuts, with zero extensions, are the shortcuts connecting the inner nodes 
to each other. 

In the absence of external force, Ro — 1. For a given F we have the largest extension in 
the bond connecting node 1 to N. Hence by increasing F the first avalanche occurs at Fi = 
kN/2. At this step we have Mi = M -l, -Ri = 2 and K x = fc(JV-2)/2. A simple calculation 
shows that for all springs connected to nodes 1 and N, we have A x = F/(k(N — 2)). So 
by increasing F we reach F 2 = k(N — 2) in which 2(iV — 3) shortcuts are teared. In this 
step, we will have a complete network consisted of all nodes between 2 and N — 1. But we 
know already that the spring connecting node 2 to N — 1 can only stand forces lower than 
k(N — 2)/2 which is much smaller than F 2 . Thus, in the second avalanche, all shortcuts will 
be removed. So F c = k(N — 2) is the threshold F in which CEED has an abrupt change. 
Next we consider the simple case of k' — 0. It is clear that in this situation, the shortcuts 
have no contribution in the elastic properties of the network and we have a trivial problem 
in this respect. But in the process of increasing F, the effect of shortcuts in quantities 
like CEED is still important and nontrivial. It is not difficult to show that in this limit 
P(A X ) oc (FN/k — A x ). Moreover, as a function of the number of avalanches a, we have 
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FIG. 3: g(y) = 2K/(kpN) vs y = pN 2 for networks of size N = 100,500 (averaged over 1000 
realizations) and N = 1000 (averaged over 300 realizations). The line represents the curve 1 + 2/y. 



R„ = l + l- A 



K„ k 



N — 1 

where [x]+ denotes the smallest integer larger than or equal to x. In Fig. El we have shown 
R versus a. Here again the relation between force and CEED is linear. If we take a and R a 
as continuous variables and use Eq. (jSJ) we obtain 

P(AR) ~ (AR)- 3 / 2 . (9) 

This behavior is also seen in mean field models of fiber bundles for the size distribution of 
avalanches Q|. 

For < k' < k, P(A X ) still consists of three delta peaks. Numerical simulations show that 
only for k! <^k the abrupt change of CEED is replaced by a continuous one (more precisely 
a staircaselike behavior). For example, if we choose N = 100 and k — 1 we obtain a nearly 
continuous transition only for k' < 0.01. We found that this value of k' is a decreasing 
function of network size. Thus for iV — > oo a continues transition occurs only at k' — 0. 



V. RANDOM SMALL- WORLD NETWORKS 



We take a one- dimensional chain of N nodes with N — 1 regular bonds. Then with 
probability p, we connect any two nodes of distance larger than 1 by a shortcut. In the 
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FIG. 4: Extension distribution for RSWNs of size N = 1000 and k = k' = 1 averaged over 5000 
realizations. The line displays a power law of exponent —2. 

following, we consider three different cases (i) k = k', (ii) k' = and (iii) < k' < k. 
A. The case k = k' 

For p > 1/N we expect to have K = kpN/2 which in the limit of p — 1 is equal to the 
stiffness of a complete network. In other word if p > 1/N, we have effectively a complete 
network in which k has been replaced with kp. On the other hand for p = and large N we 
have K = k/N. Using these limiting cases we suggest that 



Note that for large N, pN 2 /2 is equal to the average number of shortcuts. Numerical 
simulations shown in Fig. El support this scaling relation. 

As p decreases, we observe a crossover in P(A X ) from a three-modal behavior to a scale 
free one, see Fig. HJ When p > 1/N this distribution has three broad maximums instead 
of three delta peaks in a complete network. The broadening of these peaks is due to the 
random structure of the SWN. Just in the small-world regime, P(A X ) is given by a power 
law distribution of power —2. In this respect the RSWN belongs to the universality class of 
the hierarchical model introduced in section IIHI 

In our quasistatic process of increasing F, CEED changes nearly abruptly at F c . This 
behavior has been shown in Fig. EJthat displays P(AR) for single realizations of the process. 
We observe that after a few events of small size, we have a large change of order N in R which 



K kp 
N = Y 
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FIG. 5: P(AR) for single realizations with different p when k = k' = 1, N = 100. There are a few 
small changes in R followed by a large change of order N happened at critical force F c . 

we interpret it as a discontinuous behavior of CEED. Again for p 3> 1/iV, F c = kp(N — 2), 
see Fig. EJ For p — > 0, when there is at least one shortcut in the network, F c ~ kS c . These 
limiting behaviors of F c with p have been shown in Fig. El We find that average number of 
shortcuts emanating from a node, pN, determines F c . As far as the shortcuts do not overlap 
with each other, F c dose not depend on their number. Actually for a finite number of 
shortcuts, a typical configuration has almost nonoverlapping shortcuts. In this situation we 
have F c ~ k5 c . Thus we expect that F c changes considerably only when pN ~ 1. Moreover, 
from the linear feature of the system we also expect that F c should be proportional to 5 C 
and k. Thus, we suggest F c = k6 c h(pN), where h(pN) is a dimensionless scaling function. 
In Fig. El we show that this relation works well for two small-size networks and within our 
numerical errors. 



B. The case k' = 

When kl = the shortcuts have no contribution in the elastic properties of the network. 
As indicated in the study of complete networks, P(A X ) is a linear decreasing function. But 
notice that during the quasistatic process, tearing of shortcuts leads to considerable varia- 
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FIG. 6: F c vs pN in RSWNs with k = k' = 1 for N = 50 and N = 100 averaged over 5000 and 
1000 realizations respectively. The solid line displays F c = kpN. In each step of the quasistatic 
process we increase F by 0.1. The inset shows more clearly the linear behavior of F c for large pN. 
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FIG. 7: The CEED vs F in RSWNs of size N = 1000, k' = and k = 1 after 5000 (for p = 
0.1,0.001) and 10000 (for p = 0.0001) realizations. 

tions in R. In Fig. Owe show variation of R versus F obtained by numerical simulations. 
There is no linear relation between R and F. Indeed for a given F, all shortcuts connecting 
two nodes of distance larger than Lp = k/F have been already teared. This introduces an- 
other length scale. We expect R to be a function of the number of shortcuts in a subnetwork 
of size Lp, that is, 

5 = f(P L F)- (11) 
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FIG. 8: f(y) = R/N vs y = pL 2 F for RSWNs of different sizes (averaged over 10000 realizations). 
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FIG. 9: P(AR) for RSWNs of size N = 1000, k' = and k = 1. The lines display P{AR) oc 
(AR)~ a . Number of realizations are the same as those of Fig. El 

Numerical simulations presented in Fig. |H] confirm this scaling relation. 
Figure 01 displays P(AR) for various values of p. For a small value of p = 0.1, P(AR) oc 
(Ai?)~ 3 / 2 like a complete network. For smaller p, we find a power law distribution of lower 
exponent. For a very small p, P(AR) is nearly a constant function with an exponential tail. 

C. The case < k' < k 

As before we take k — 1 and decrease k! from one. Numerical simulations show that 
in this case the network stiffness is proportional to k! . The effects of k! and p on elastic 
properties are the same. It means that with respect to elastic properties, the decreasing 
of k! at a fixed p is equivalent to the decreasing of p at a fixed k! . In Fig. El we display 
the extension distribution when k' = 0.1. It is observed that when p > 1/N we have a 
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FIG. 10: Extension distribution for RSWNs of size N = 1000, k = 1 and k! = 0.1 averaged over 
5000 realizations. 
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FIG. 11: Correlation coefficient of Ax and b for N = 100 after averaging over 5000 realizations. 

multimodal distribution in contrast to the three-modal distribution of case k = k' for the 
same value of p. Recall that we did not have such a behavior for lower values of p when 
k = k' . Indeed this new multimodal behavior is observed only when we have a large number 
of shortcuts. 

The extension distribution of RSWNs shows that there are some shortcuts bearing very 
large forces compared with the other ones. Certainly this is because of their essential roles 
in the network structure. A good measure of centrality of bonds in a network is their 
betweenness |2j- Suppose that we have n shortest pathes connecting node 1 to N. A given 
shortcut may contribute in n s of these pathes. Then, the betweenness of this shortcut is 
defined as b = n s /n. We define 



< bA x > 



< b >< A x > 



TbA. '■= 7== > ( 12 ) 

as a measure of correlation between extension and betweenness of shortcuts. In this def- 
inition a^ x and cr b are variances of A x and b, respectively. In FigJTT] we show how the 
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FIG. 12: P(AR) for a network with N = 1000, p = 0.001 and k = 1. The data are results of 
averaging over 10000 = 0) and 5000 (A/ = 0.1,0.2) realizations. The lines represent power laws 
of exponent —1 and —3/2. 
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FIG. 13: F max vs k' for TV" = 100 and k = 1 (averaged over 1000 realizations). The inset shows more 
clearly the linear behavior of data for large k! . The lines display linear functions F max = ao + a\k' . 

correlation coefficient depends on k! and p. As expected, r b ^ x has a considerable positive 
value in the small-world regime. The correlation coefficient is nearly independent of k' ex- 
cept for a rapid decrease to zero for k' — ► 0. 

In the quasistatic process, as it happens in complete networks, a considerable number of 
avalanches takes place only for a significantly small k! . Indeed, the larger p the lower value 
of k'\ we need to see a continuous behavior of R. In Fig. ^]we compare P(AR) for two 
cases of k! = and k' = 0.1. The figure shows that as expected, by increasing k', P(AR) ap- 
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proaches to that of a complete network. Indeed, in the small-world regime and for a nonzero 
k', P(AR) obeys a power law of exponent —a with a < 3/2. This exponent approaches to 
—3/2 as we increase k' or p. 

In this case it is not easy to define k' c (p), the value which separates continuous and discon- 
tinuous regimes. Instead we calculate F max , the force which leads to R max — N — 1. We 
expect that a change in the behavior of F max with respect to k' signals a crossover from 
a continuous to a discontinuous region. In Fig. EH we show the variation of F max versus 
k' obtained by numerical simulation. We see that by decreasing k', the linear behavior of 
Fmax changes and then it saturates for k' — > 0. This crossover occurs for a lower k' as one 
enhances the number of shortcuts. 

VI. CONCLUSIONS 

We studied the static and quasistatic properties of SWNs. We observed that the 
network's structure significantly affects the static and quasistatic behavior. Thus, we 
could draw some conclusions about the structure of networks through the study of these 
properties. 

The summary of the main results for RSWNs have been represented in table |l] and Fig. 
[T41 We found that for very small values of p, the effective stiffness of the network is 
comparable with that of a complete network. When k = k' the extension distribution of 
shortcuts is a power law of exponent —2. In this respect RSWNs behave like a hierarchical 
network introduced in this paper. There was also a strong positive correlation between 
the betweenness of a shortcut and its extension. It means that just by looking at the 
distribution of extensions in a network, one could be able to distinguish which shortcuts 
are more central. 

In the quasistatic part, we showed that by increasing F, CEED could have a continuous or 
discontinuous transition. In general, to have a continuous transition we need a much smaller 
spring constant for shortcuts rather than regular bonds. In the case of a discontinuous 
transition, the critical force is determined by the average number of shortcuts per node. It 
was found that for a continuous transition, P(AR) is given by power law distributions of 
the exponent —a with a < 3/2. 
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FIG. 14: Schematic representation of phase diagram. Dashed region indicates where we have a 
continuous behavior of R vs. F. For p < 1/N 2 (black region) we can not speak of continuous or 
discontinuous behavior of R. Note that as N approaches to infinity the continuous region reduces 
to the line k' = 0. The inset shows the variation of a with p when N = 1000, k = 1 and k' = 0. 
The errors in a are of order 0.01. The same behavior will be observed for a when we fix p and 
increase k'. The solid line in inset indicates the expected value of a for large p i.e. 3/2. 





P(A X ) 


P(AR) 


k = k' 




5ar,n (an abrupt change of R) 


k! = 


a - aiA^ 


(AR)- a (a < 3/2) 


< k' < k 


(Aj;) _/3 (/3 < 2)+ oscillations (for large p) 


1 


(AR)~ a (a < 3/2), k' < k' c (p) 
5ar,n, k' > k' c (p) 


\ 



TABLE I: Summary of results for P(A X ) and P(AR) in different cases. Here k' c (p) gives the line 
in (p, k') plane that separates continuous and discontinues behavior of CEED. 

Acknowledgments 

We would like to thank S. N. Rasuli for useful discussions. We are grateful to MirFaez 
Miri for comments and careful reading of the manuscript. 



15 



[1] R. Albert and A. -L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). 

[2] M. E. J. Newman, SIAM Review 45, 167 (2003). 

[3] D. J. Watts and S. H. Strogatz , Nature 393, 440 (1998). 

[4] L. A. N. Amaral, A. Scala, M. Barthelemy, and H. E. Stanly, Proc. Natl. Acad. Sci USA 97, 
11149(2000). 

[5] A. Barrat, M. Weigt, Europ. Phys. J. B 13, 547 (2000). 

[6] L. de Arcangelis, H.J. Herrmann, Physica A 308(1-4), 545 (2002). 

[7] A. D. Sanchez, J. M. Lopez, M. A. Rodriguez, Phys. Rev. Lett. 88, 048701 (2002). 

[8] C. Castellano, D. Vilone, A. Vespignani, Europhys. Lett. 63, 153 (2003). 

[9] E. Almaas, R. V. Kulkarni, D.Stroud, Phys. Rev. E 68, 056105 (2003). 

[10] K. Medvedyeva, P. Holme, P. Minnhagen, B. J. Kim, Phys. Rev. E 67, 036118 (2003). 

[11] L. G. Morelli, G. Abramson, M. N. Kuperman, Eur. Phys. J. B 38, 495 (2004). 

[12] S. Jespersen, I. M. Sokolov and A. Blumen, Phys. Rev. E 62, 4405 (2000). 

[13] S. Jespersen, I. M. Sokolov, A. Blumen, Arxiv preprint |cond- mat / 0004392 . 

[14] S. Jespersen, A. Blumen, Phys. Rev. E, 62, 6270 (2000) . 

[15] M. Vendruscolo, N. V. Dokholyan, E. Pad, M. Karplus, Phys. Rev. E 65, 061910 (2002). 

[16] Y. Kantor, M. Kardar, Phys. Rev. E 54, 5263 (1996). 

[17] P. Sen, B. K. Chakrabarti, J. Phys. A: Math. Gen, 34 (38), 7749 (2001). 

[18] D. Sornette, J. Phys. A 22, L243 (1989). 

[19] PC. Hemmer and A. Hansen, J. Appl. Mech. 59, 909 (1992). 



16 



